Interference patterns of Bose-condensed gases in a two-dimensional optical lattice 
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For the Bose-condensed gas confined in a magnetic trap and in a two-dimensional optical lattice, 
the non-uniform distribution of atoms in different lattice sites is considered based on Gross-Pitaevskii 
equation. A propagator method is used to investigate the time evolution of 2D interference patterns 
after (i)only the optical lattice is smithed off, and (ii)both the optical lattice and the magnetic trap 
are smithed off. An analytical description on the motion of side peaks in the interference patterns is 
presented by using the density distribution in a momentum space. 
PACS number (s): 03.75. Ft, 05.30. Jp 

I. INTRODUCTION 

In the last few years, due to the experimental realization of Bose-Einstein condensation in 1995 [1], we witnessed the 
remarkable experimental and theoretical advances in the study of ultra-cold Bose gases [2] . Recently, a Bose-condensed 
gase confined in an optical lattice (created by retroreflected laser beams) and in a magnetic trap is investigated (see 
[3] and references therein). In such system, many subcondensates form and they are confined in the minimums of 
optical lattice potential. This provides an ideal system for investigating many interesting properties of Bose-Einsten 
condensate (BEC), such as the quantum transition from a superfluid to a Mott insulating state of a ultra-cold Bose 
gas [4]. 

The coherence property of a BEC confined in an optical lattice and in a magnetic trap can be obtained by observing 
the time evolution of the interference patterns when both the optical lattice and the magnetic trap are switched off. 
The interference patterns obtained in this way is similar to the diffraction patterns of a coherent light from a grating. 
For the BECs in one-dimensional (ID) [5], 2D [6] and 3D [4] optical lattices, the interference patterns are clearly 
shown experimentally. In [4], the quantum transition from a superfluid to a Mott insulating state is confirmed by 
observing the characters of the interference patterns. 

In a recent experiment, Miiller et al [7] studied the time evolution of the interference patterns when only a ID 
optical lattice is switched off. It is shown that, due to the presence of the magnetic trap, there is a periodic harmonic 
motion in side peaks. This gives us a new opportunity to investigate the properties of coherent matter wave such as 
the collision between the side peaks. In our previous work [8], a propagator method is used to study the time evolution 
of the interference patterns formed in a ID optical lattice, and the harmonic motion of the side peaks predicted is in 
agreement with the experimental result given in [7]. In the present work, we generalize the method developed in [8] 
to a 2D case, i. e. we consider the time evolution of the interference patterns for a Bose-condensed gas confined in a 
2D optical lattice and in a magnetic trap, with the emphasis especially on the situation when only the optical lattice 
is switched off. 



II. A BOSE-CONDENSED GAS IN COMBINED TRAPPING POTENTIALS 
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For the Bose-condensed gas confined in a 2D optical lattice and a magnetic trap with an axial symmetry along 
z-axis, the dynamics of the condensate at low temperature is described by the well-known Gross-Pitaevskii equation 
[3] 
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where tu± and tu z are respectively the harmonic angular frequencies of the magnetic trap in the radial and axial 
directions, g — ATT 2 fi 2 a s /m with a s the s-wave scattering length. When the polarization vectors of two perpendicular 
optical standing waves are orthogonal, the 2D optical lattice potential takes the form [6] 



Vlatt — Uq 
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(2) 
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where Uq denotes the depth of the optical lattice which can be increased by increasing the intensity of retroreflected 
laser beams, while A is the wavelength of the retroreflected laser beams. We consider here the case of u)± « uj z and 
hence the condensate is disk-shaped with a larger extension in the x and the y directions when in the presence of only 
the magnetic trap. After the optical lattice potential Vi a tt is switched on a lot of subcondcnsates appear which are 
confined in the minimums (lattice sites) of the optical lattice potential. For the optical lattice created by the laser 
beams with wavelength A, d = A/2 is the period of the optical lattice potential, which can be regarded as the distance 
between two neighboring lattice sites. 

We assume that the width of subcondcnsates in each lattice site, a, is much less d, which can be realized in 
the present experiment by simply increasing the depth of the optical lattice. We consider here the case that the 
subcondcnsates in different lattice sites are fully coherent, i.e., the chemical potential of these subcondensates are 
identical and hence the condensed state wave function (or called order parameter) \1/ of the overall subcondensates 
can be written as: 

*= E ^k v (x,y,z)e- i ^ n , (3) 

iky 

where \i is the chemical potential, kj (j = x, y) denotes the fcj-th lattice site in the j-direction. 
Substituting (3) into Eq. (1), we have 
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In obtaining Eq. (4) we have used the condition that the overlap between neighboring subcondensates can be neglected 
in the case of a « d. To illustrate the character of $k x k v clearly, using the coordinate transformation x — k x d — > x, 
y — k y d — > y and making a Taylor expansion of the optical lattice potential to the quadratic terms with respect to the 
lattice site {k x ,k y }, it is straightforward to get the following equation: 

2~m ^ 2 + \ muj2 ^ x + kxd ^ + ( y + fc f rf ) 2 ) + \ muj2 z z2 + \ m ^ 2 i- ( x2 + V 2 ) 



+9 \®k x k y 



®k x k y = ^k^ky, 



(5) 



Equation (5) determines the wavefunction $fe x fe y at the lattice site (k Xl k y ). In the present experiment, the intensity 
of the retroreflected laser can be increased so that one can have ui±_ > > uj± . For the subcondensate confined in the 
lattice sites, one has — d/2 < x < d/2 and — d/2 < y < d/2. Thus we obtain x « k x d and y « k y d except for 
the subcondensate in the central lattice site. As a result the term ^moj 2 [((x + k x d) + (y + k y d) ) in Eq. (5) can be 
approximated as \mLo\{k 2 x d 2 + k 2 d 2 ). As for the subcondensate in the central lattice site, because uj± >> uj±, the 
above approximation is valid too. Thus Eq. (5) can be reduced to 



~2~^i y2 + o muj2z2 + o muj2 t- ( x2 +y 2 )+9 |*fc x fc H | 
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where ^ klc k y = M~ 
on the lattice site {k x ,k y }. Assuming that fi 



)\d 2 {k 2 + k 2 ), which can be regarded as an effective chemical potential of the subcondensate 



\mJ 2 1 _d 2 k 2 M , the effective chemical potential then reads 



Mfc x fe y = -muj\d 2 {k 2 M 
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(7) 



Because the occupied lattice sites lie within a circle, ^k 2 M can be regarded as the number of subcondensates induced 
by the optical lattice. 

In the present work, we are interested in the case Tiuj z « fi « huj±. In this situation, each subcondensate in the 
optical lattice is a cigar-shaped one along the z-direction thus it displays a quasi-lD character. So a Thomas-Fermi 
approximation can be applied in the z direction. After integrating out the fast- varying variables x and y in Eq. (6), 
one obtains [9] 



ti z d 2 



1 



2^ + o TO ^ +9iD\<Pk.k v {z)\ 



¥>k X ky (Z) = Mfe x fe y ^fe x fe y (Z) . 



(8) 



2 



where ifk x k y (z) is the wave function of the subcondensate at the lattice site (k x ,k y ) in the z-direction and gio = 
2h 2 a s /mT\ (l± = (h/mu ± ) 1/2 ) is an effective coupling constant in one dimension. 

In the z direction, using the Thomas- Fermi approximation we obtain the atomic number N kxky in the lattice site 

^y} &S 

Thus the ratio between the atomic number in the lattice site (k x , k y ) and the atomic number in the central lattice 
site (0,0) reads 

N kxky = N 00 (l - ■ (10) 

We see that the uniform chemical potential of the fully coherent subcondensates and the presence of the magnetic 
trap result in a non-uniform distribution of the atomic number in different lattice sites. It is worth pointing out that 
the atoms can move from one lattice site to another ones when the subcondensates are fully coherent (i.e. the system 
is in a superfluid state). Thus the atomic number in lattice sites discussed here should be regarded as an average one. 
In addition, the value of &m can be obtained by considering the condition N = J2k k ^k x k y - Here N is the total 
atomic number in all subcondensates, the prime in the sum means that the summation about k x , k y should satisfy 
the condition k 2 + k 2 < k 2 M . Based on this analysis, ku is given by 

k M = -. t-tt • (11) 

In the present work, the parameters used are ui± = 24 x 2n Hz, u> z = 220 x 2n Hz and uj± = 10 4 x 2n Hz 
for 7V = 10 5 87 Rb atoms. Note that uj± here is about one half of the experimental value used in [6], thus the 
subcondensates discussed here are fully coherent. It is obvious that the condition huj z « [i « hu>± used above 
holds for these parameters. 



III. THE BOSE-CONDENSED GAS IN A MOMENTUM SPACE 



We now investigate the density distribution of the Bose-condensed gas in a momentum space. Note that for the 
parameters used here, the Thomas-Fermi approximation is no more valid in the x and y directions. For a larger value 
of J7o it is reasonable to assume that the density distribution of subcondensate in the optical lattice has a form of 
Gaussian distribution in x-y coordinate space, i.e. $o (x,y) ~ ^2 k k exp[— ((x — k x d) 2 — (y — k y d) 2 ) /2a 2 }. Due to 
the periodicity of the lattice, the wave function in momentum space takes the form 

sin[(2fc M + l)p x d/h] sin[(2fc M + l)p y d/h] 
smp x d/2n smp y d/2n 

where $o(p x ,p y ) ~ exp[— (p 2 + p 2 )a 2 /2h 2 ]. We see that the momentum distribution of the Bose-condensed gas 
exhibits sharp peaks at pj = 2wrijh/d (j = x,y). The width of each peak in momentum space can be approximated 
as Apj ~ irh/ (2&m + 1) d. From this character one can obtain directly the following results: 

(i) As pointed out in [5,8], the density distribution of the subcondensates in the momentum space provides important 
information for the evolution of interference patterns. The sharp peaks of the momentum distribution of the condensate 
with nonzero rij implies that there would be side peaks in coordinate space after the optical lattice is switched off. 
The classical approximation of the initial velocity Vj — Pj/m can give a quite well description for the motion of the 
side peaks. When only the optical lattice is switched off, due to the presence of the magnetic trap the harmonic 
motion of the center of mass of the side peaks is given by 

2im x Ti . , . 2irn v Ti . , . , . 

r nx ,n y = -1 sin (ujj_t) e x H a — sin (u±t) e y , (13) 

where e x and e y are two unit vectors along the x and y directions. When both the optical lattices and magnetic trap 
are switched off, the motion of the side peaks is given by 
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2imxht 2Tm y ht 

Tn *' n » = ~^d~ Gx + ^d~ ey ' (14) 

(ii) From the density distribution in momentum space n (p x ,p y ) — |$o (Px,Py)\ 2 i the relative population of the side 
peaks with respect to the central peak can be approximated by exp[— 47r 2 (n 2 + n 2 )a 2 /d 2 ]. This formula is very useful 
when experimental parameters are chosen to observe clearly the side peaks in the interference patterns. 

IV. DENSITY DISTRIBUTION OF THE CONDENSATE AFTER SWITCHING OFF ONLY THE 

OPTICAL LATTICE 

In [8], evolution of interference patterns of the subcondensates induced by a ID optical lattice is investigated when 
only the optical lattice is switched off. In this section, we generalize the method developed in [8] to consider the 
evolution of interference pattern when swithching off only the 2D optical lattice. Using the density distribution in 
coordinate space and the result given by Eq. (10), the normalized wave function of the Bose-condensed gas at time 
t = takes the form 

(x, y, Z ,t = 0) = A n ^k y (*> V,t = 0) ®k x k y (z, t = 0) , (15) 

where A n = yJb/2/irakM is a normalization constant. In the above equation, $k x k y (z, t = 0) represents the normal- 
ized wave function in the z— direction, and ®h x k v {x, y,t = 0) is given by 

/ u2 , 1.2 \ 3/4 

tf,t = 0) = h-^pj e -((*-W) 2 +(,-M)W. (is) 

Note that the factor (l — {k 2 + k 2 ) /k 2 J f'^ in the above expression reflects the non- uniform distribution of the atoms 
in the lattice which is given by Eq. (10). 

When the optical lattice is switched off, the evolution of a interference pattern can be illustrated by considering 
the density distribution n(x,y,t) in x-y coordinate space: 

n(x,y,t) = N J \^{x,y,z,t)\ 2 dz = N \V (x,y,t)\ 2 . (17) 

In the case of the effective chemical potential [ik x k y being much smaller than the ground state energy huj±/2 of the 
atoms in the lattice site (k x ,k y ), the noninteracting model can provides a good description [5,8]. When neglecting the 

interatomic interaction, the normalized wave function ^ (x, y, t) = A n J2k x k y ^k^k y {x, y, t) can be obtained through 
the well-known propagator method [8]. Recently, the vortex dynamics of parabolically confined Bose gases is also 
investigated by the propagator method [10]. Using this technique the condensed state wave function at time t can be 
obtained by using the following integral equation [11]: 

/oo 
K(x,y,t;x!,yi,t = 0)* (xi,yi,t = 0) dxidyi, (18) 
-oo 

where the propagator K(x, y,t;x\,yi,t — 0) is given by 

K(x,y,t;x 1 ,y 1 ,t = 0)= JJ K^Tj^rjut = 0). (19) 

In the above equation, r x = x (r y = y), r x \ = X\ (r y \ = y{), and 

1/2 



K j (r j ,t;r jl ,t = 0) = 



muj± 



6XP {2^^ ^ + ^ ^ - ^ } ' (20) 



_2tt ih sin cu±t_ 

By Eqs. (16) and (18) it is straightforward to get the following analytical result of ^ (x, y, t): 
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(x, y,t) = A n > l 



where 



sinwjj (cotw_i_i + 27) 



fr 2 



exp 



n s j ' *) • 



(21) 



3=x,y 



(kjdcosLo^t — rj) 



exp 



i(kjd cos u>i_ t — rj) cotu^ i 
27CT 2 sin 2 wjj (cot 2 w_i_i + 7 2 ) 



exp 



2er 2 sin 2 ujxt (cot 2 w ± t + 7 2 ) 

i (r 2 cosujj + fc 2 d 2 cosw^t — 2kjrjd) 
2-fa 2 smu±t 



(22) 



where the dimensionless parameter 7 = h/muj±a 2 . 

From the result given by Eq. (21), we see that the period of \& (x, y, t) is T = 2tt/uj±, while the period of the density 
distribution n(x, y, t) is it/uj±. The density at x = 0, y = reaches a maximum value at time t m = (2m — 1) tt/ujj_ 
with m a positive integer. At time t m , |\& (x = 0,y = 0, t m )\ is given by 



\y(x = 0,y = 0,t v 



1 2 2 

I =OL- 



4 2 

7 2 



47rfc M 
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(23) 



According to the fact that there is only a small number of atoms in each lattice site, we assume here that the atoms 
in a subcondensate are in the ground state of the corresponding lattice site. Thus we have a — y / h/2m,LUj_ and it is 
easy to verify that for the parameters used here satisfy a « d. This shows that in the derivation of Eq. (4) the 
condition a « d holds. From the result given by Eqs. (21) and (22), one can obtain the density distribution of the 
interference pattern. Displayed in figure la is the density distribution n(x,y,t) (in unit of -/VqA 2 ) at t = 0.3it/uj±. 
Note that in all figures plotted in this paper, the coordinates x and y are measured in the unit d. The central and side 
peaks are clearly shown in the figure. Due to the presence of the magnetic trap, there is a harmonic motion of the 
side peaks. Displayed in figures lb-d are the central and side peaks of the interference pattern at t = 0.5tt/lj±. From 
the result shown in figure lb, we see that the central peak is a very sharp one, analogously to the case of ID optical 
lattice [8]. Shown in figure 2 is the evolution of the interference pattern at different times. The evolution in the width 
of the central and side peaks and the motion of the side peaks is obvious. Due to the presence of the magnetic trap, 
there is no obvious evolution in the z-dircction. 

V. DENSITY DISTRIBUTION OF THE CONDENSATE AFTER SWITCHING OFF THE COMBINED 

POTENTIALS 



Wc now turn to discuss the time evolution of the interference patterns when both the magnetic trap and the optical 
lattice are switched off. In this case the propagator can be obtained by setting ui± — > in Eq. (19). Thus we get 



K b (x,y,t;x 1 ,y 1 ,t = 0)= JJ_ K hj {rj,t;rji,t = 0), 



where 



K bj (rj,t;rji,t = 0) 



The analytical result of $ (x, y, t) is then 
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=x,y 


2iriht 


-1/2 j 




exp I 


m 





im (rj — TjiY 
2M 



(24) 



(25) 



V(x,y,t)=A n J2 1 - 



k 2 



(26) 



where 
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^to'^U+ieJ exp 

where the dimensionless parameter = Tit/ma 2 . Shown in figure 3a is the density distribution (in unit of N A n ) at 
t = 0.3ir/w±. The density distribution at t = 0.57r/cjj_ is given in figure 3b. We see that the peaks are wider than 
those in the case when only the optical lattice is switched off. Displayed in figure 4 is the evolution of the interference 
pattern with the development of time. 

As for the evolution in the z-direction, analogously to the case of the ID optical lattice [5] one can obtain that the 
width R z (t) of the central peak is described by the asympotic law R z (t) = R z (0) uj z t after the central peak is formed. 



(rj - kjd) 
2a 2 (1 + 6 2 ) 



exp 



i (rj — kjd) 
2a 2 (0 + 1/6) 



(27) 



VI. DISCUSSION AND CONCLUSION 



The time evolution of the interference patterns of a Bose-condensed gas confined in a 2D optical lattice and a 
magnetic trap has been investigated by using a propagator method. Based on this method the analytical results of 
the condensed state wave function and motion of the side peaks have been provided. When the effective chemical 
potential in each lattice site is much smaller than Hoj±, a noninteracting model can give a quite well description 
for the evolution of the interference patterns. Nevertheless, the interaction between atoms can give an important 
correction to the central peak when only the optical lattice is switched off. After switching off the optical lattice, 
at time t m the central peak is very sharp and hence the interaction between atoms can not be neglected. Generally 
speaking, starting directly from a general mean field theory, one can give the interaction correction to the density 
distribution of the central peak. However, this would be a challenging work due to the fact that there is not a single 
phase for the central peak due to the evolution of the wave packet. In [8], for a ID optical lattices a simple method is 
developed to investigate the interaction correction by using a Gaussian density distribution of the central peak and 
energy conservation. In the case of the 2D optical lattice, the analysis is similar and the interaction correction to the 
maximum density of the central peak at t m can be calculated straightforwardly. Assuming that (3 denotes the ratio 
of the density at x — 0, y = between the weakly interacting and noninteracting interference patterns, at t m we have 



P= T- p 7 wT^, (28) 

1 + bint {a±-ideal) / &all 



where a±-id ea i is given in Eq. (23), with 



h 2 1 

Eall = N [ + -mti\o 2 ) , (29) 



and 



bint (Oi±-ideal) = ^T, : • (oUJ 

For the parameters used here, a simple calculation shows that (3 — 0.43. We see that quite differently from the case 
of ID optical lattices (see Rcf. [8]), the weak interaction between atoms results in a change in a deep way for the 
density of the central peak. This is not a surprising result due to the fact that in the case of the 2D optical lattice, the 
central peak is cigar-shaped, while the central peak is disk-shaped for ID optical lattices. Another possible interaction 
correction to the central peak at time t m would be a four-wave mixing which deserves a detailed study in future. 

In addition, the interaction between atoms plays an important role in the colliding process of the side peaks, and 
this would be an interesting problem deserving to be explored further. In the case of the Bose-condensed gas in a ID 
optical lattice, the collisional haloes have been observed in the experiment [7]. When a 2D optical lattice is used to 
confine an axial symmetric Bose-condensed gas, there would be a collision between four side peaks, rather than the 
collision between two side peaks as in the case of ID optical lattice. This gives us new opportunity to investigate 
the collision between side peaks both theoretically and experimentally. The method developed here can also be used 
straightforwardly to study the cigar-shaped Bose-condensed gas in a 2D optical lattice. For the cigar-shaped Bose- 
condensed gas, the harmonic angular frequencies of the magnetic trap is u) x — uj z » u y . When the 2D optical lattice 
is imposed in x and y directions, after the optical lattice is switched off the motion of the side peaks is given by 



2irn x Ti . . , 2im v Ti 

sin (w x t) e x H *— 

y muj±d muj±a 



r n x ,n„ = ; sin (w x t) e x + — v " sin (uj y t) e y . (31) 
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We see that the motion of the side peaks in x and y directions is quite different because lu x >> u> y . fn this case, 
the side peaks in the x-direction has a collision before that in the y-direction, and this gives us an opportunity for 
observing the collisions one after another. It is obvious that the method developed here can also be used to investigate 
the interference patterns of the Bose-condensed gases in a 3D optical lattice which has been considered recently in 
[12] by using a different method. 
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Figure 1: The interference patterns after only the 2D optical lattice is switched off. Figure la shows the density 
distribution n(x,y,t) (in the unit N A^) at t = 0.3tt/lj±. Figures lb-d show the central (figure lb) and side peaks 
(figure lc for {n x = 1, n y = 0} and figure Id for {n x = 1, n y = 1}) at t = 0.5it/uj±. In all figures, the coordinates x 
and y are measured in unit of d. Due to the presence of the magnetic trap, the central and side peaks are very sharp 
at t = 0.5tt/uj±. 

Figure 2: The evolution process of the central and side peaks after only the 2D optical lattice is switched off. The 
interference patterns are shown for different time-of-flight t = 0, 0.1tt/luj_, 0.2ir/ui±, 0.3ir/ui±, 0Air/ui±, and 0.5it/uj±. 

Figure 3: The interference patterns after both the 2D optical lattice and the magnetic trap are switched off. Figures 
3a and 3b show the interference patterns for t = 0.3tt/lli± and 0.5tt/uj±, respectively. The density distribution n (x, y, t) 
is measured by NqA^, while the coordinates x and y are in the unit of d. 

Figure 4: The evolution process of the central and side peaks for for different time-of-flight t = 0, 0.i7r/wj_, 0.27r/a>^, 
0.37r/w^, 0Air/oj±, and 0.5tt/lj± after both the 2D optical lattice and the magnetic trap are switched off. 
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